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PERTURBATIONS OF THE INTERPLANETARY MAGNETIC 
FIELD BY THE LUNAR WAKE 


Abstract 

Detailed measurements of the interplanetary magnetic field in the 
vicinity of the moon have been performed on the Explorer 35 selenocentr ic 
spacecraft launched July 1967. The observational evidence fails to 
reveal the existence of a lunar bow shock wave. A pseudo-magnetosphere 
as suggested by Gold and interpreted from Luna 10 measurements is not 
observed. The interplanetary magnetic field appears to be convected 
through the lunar body without a large scale distortion of its direction 
or magnitude. Perturbations as much as 307 o of the average magnitude 
are noted in the solar plasma umbra and penumbra. As the satellite 
passes through the leeward flow an alternating pattern of magnitude 
increases and decreases is observed in the penumbra while generally 
only an increase is observed in the umbra. Using a theoretical model 
of plasma flow due to Whang, a first order solution of the perturbed 
interplanetary magnetic field is compared with observations. It is 
concluded that the perturbations can be partially explained on the basis 
of the magnetization, gr’ad-ient and curvature currents induced in the disturbed 
solar plasma flow. The umbral increase and the innermost penumbral 
decrease are consistent with the first order theory and it is suggested 
that a higher order approximation is required to explain the newly 
detected penumbral increases and additional penumbral fluctuations. 
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Introduction 

A previous report (Ness et al., 1967) has presented early results 
from the NASA-GSFC magnetic field experiment carried on Lunar Explorer 
35. This satellite was placed into selenocentr ic orbit on July 22, 1967 
with aposelene = 9388 + 100 km (5.4 R^) , periselene = 2568 ± 15 km 
(1.5 R^) , period=11.5 hours and initial aposelene-moon-sun angle = 304°. 
The initial interpretation of these data indicated the absence of a 
bow shock wave formed in the supersonic flow of solar plasma past the 
moon. The NASA-GSFC magnetic field experiment does not yet find a bow 
shockwave across which the interplanetary magnetic field increases 
abruptly by a factor of 2 to 4 nor the development of a fluctuating 
magnetic field regime in a sheath layer. In a related study, Lyon et 
al. (1967) reported the existence of a plasma shadow region behind the 
moon in which a detectable solar plasma flow was absent. In addition 
there appeared to be no evidence in these data for a shock wave which 
would be detectable by the presence of a sheath layer of "hot" electrons 
and a disturbed proton flux. These results do not support the earlier 
interpretations of the existence of a lunar magnetosphere from measure- 
ments performed by the Luna 10 spacecraft in 1966 (Dolginov et al . , 1967). 

A feature of the magnetic field observations noted by Ness et al. 
(1967) and Colburn et al. (1967) was the existence of increases in the 
interplanetary magnetic field magnitude in the solar wind shadow region 
and associated decreases of the magnitude in the penumbral region of 
the solar plasma flow. On the basis of the absence of a large scale 
disturbance of the interplanetary magnetic field in the flow past the 
moon. Whang (1967) developed a theory of plasma flow based upon a guiding 
center approximation. In this model the solar plasma is treated as a 
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collisionless gas with one dimensional thermal motion of the 
protons only along magnetic field lines. A long non-axially 
symmetrical wake region is developed behind the moon in which a plane 
of symmetry exists defined by the local solar wind velocity and inter- 
planetary magnetic field direction. The exact geometry of the wake 
depends upon both the direction of the interplanetary magnetic field 
and the thermal energy of the protons, i.e., the ratio of the mean 
random velocity of the solar plasma to its bulk velocity. 

It is the purpose of this paper to present additional observational 
evidence supporting the original conclusions by Ness et al. (1967) 
on the absence of a detectable shock wave and to study in detail the 
perturbations of the interplanetary magnetic field in the lunar wake. 
These measurements will then be compared with theoretical results 
based upon the Whang model of the solar plasma flow. The perturbations 
are due to magnetization, gradient and curvature currents 
in the disturbed solar plasma flow. Studies of the solar plasma flow 
around the moon have been conducted by Michel (1967 a,b) and Johnson 
and Midgley (1967) and their relation to the present study shall be 
discussed. 

The data to be presented in this report cover only the interval 
from July 29“ August 6, 1967. These data follow shortly after those 
discussed in the earlier publication by Ness et al. (1967) and represent 
an interval of time during which the moon*s position relative to 
the earth changes from third quarter to new moon. During this time 
interval RF shadowing of the spacecraft is restricted to near the time 
of aposelene and a continuous monitoring of the interplanetary magnetic 
field in the leeward portion of the lunar wake is possible. Analysis 
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of data obtained from subsequent orbits leads to the same conclusions 
reached from a detailed study of this limited interval of data. As 
will be seen there exists a variety of new observational features, 
most of which fit into a coherent pattern of small perturbations. 



Data Analysis 
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Details of the NASA-GSFC magnetic field experiment on Explorer 35 
have been presented in the previous publication by Ness et al. (1967). 

A particularly interesting but also complicating aspect of the spacecraft 
system has been its performance in the optical shadow of the moon. In 
the absence of a solar input, all directional references on the 
spacecraft would be lost were it not for the inclusion of a pseudo-sun 
pulse which provides a timing reference signal through the lunar shadow. 

The time interval between successive sun pulses during the lunar shadow 
is equal to the last spin period obtained in solar 

illumination. If the spacecraft rotation rate remained constant through 
lunar shadow then this would provide a reasonable estimate for 
computation of directional references. 

An intrinsic feature of the spacecraft, however, is that in lunar 
shadow the rotation rate increases, due to thermal contraction and the 
conservation of angular momentum. As the spacecraft loses the heat 
input from the sun it is cooled and contracts. This contraction 
decreases the angular moment of inertia and hence the spacecraft rotation 
rate must increase. As the spacecraft emerges from lunar shadow the 
process is reversed and the spacecraft de-spins to its original 
rotation rate. 

The net change in rotation period throughout this optical shadow, 
of generally 30 to 40 minutes in length, is only (a few) about one part in a 
thousand, actually 2 parts in 1837. However, during the 30 minute 
interval, with a nominal 2.29 second spin period, the slight change is 
sufficient so that directional measurements at the end of the optical 
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shadow are in error by 360° or more due to the cumulative effects 

during the more than 1000 rotations of the satellite. 

A detailed study of this phenomenon and a numerical-analytical 

technique for rectification of directional measurements in the optical 

shadow of the moon has been conducted by Taylor (1967). In this paper 

experimental results are presented in which directional references 

have been corrected for the variable spin rate of the spacecraft both 

in optical shadow and thereafter. The accuracy of the correction 

o 

procedure is estimated at + 10 . It may be noted that this correction 

is applied only to the azimuthal angle, 0, so that none of the directional 

uncertainty is present in the latitudinal angle, 0. 

While in the direct solar illumination directional measurements 

are limited by the sensitivity and accuracy of the instrument (+ .ly 

and ± 0.3y, respectively) and the accuracy with which the spin axis 

orientation itself is known (+ 2°). In a typical interplanetary magnetic 

field of 6 gamma the sensitivity uncertainty leads to a directional 

o 

uncertainty of + 1 . Thus for all the measurements reported here, 
magnitudes are presented with an accuracy of + 0.3y and directions are 
presented with an uncertainty of less than 3° while the spacecraft is 
in sunlight, but perhaps as much as 10° during the optical shadow. 

Experimental data to be reported are the sequence averages in 
selenocentric solar ecliptic coordinates computed from nominally 16 
measurements at 5.12 second intervals. Computation of the RMS 
deviation of both the magnitude (6F) and the separate components (6C) 
has been included. The magnitude and component deviations are defined 


as 
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, N „ i 

6F = [— i- E (F l - <F>) 2 ] (1) 

N i=l 

and 

6C = S [(B xi - <B x >) 2 +(B yi - <By>) 2 +(B zi - <B,>) 2 ] 

N < 16. These two quantities are particularly diagnostic in a 

study of the solar wind interaction with the moon for they are sensitive 

to the existence of a sheath region in which fluctuating magnetic 

fields would be readily identified. Note that 6C > 6F always. 

The relative geometry of the satellite orbit has been shown in 

Figure 2 of the earlier publication of Ness et al . (1967). The data 

to be presented in this paper cover the interval from July 29 -August 6 

and the relative orbital parameters are changed only slightly in a 

selenocentric solar ecliptic system, the aposelene-moon-sun angle 

decreasing from 0ggg = 290° to 284°. Since the orbital plane of the 

satellite is inclined at 169° to the selenocentric solar ecliptic 

equatorial plane, i.e., almost parallel, the three dimensional 

positional reference can reasonably be approximated by the selenocentric 

solar ecliptic longitude (0 ) and the selenocentric distance RAD 

SSE 

as measured in units of lunar radii (1 = 1738.1 km). For 

interpretation of the results, the position of the orbit relative to 
the previously defined "plane of symmetry" is also important. 

The data to be presented are final production data from which possible 
noise points have been eliminated. The threshold criteria used to delete 
a particular average field vector are 

( i) 6C/F > .33 and 6C > .93 or 

(ii) more than 207, of the individual measurements used in 
computing the average are missing, i.e. N < 12. 

There are also obvious gaps in the data plots which correspond to 


periods when the spacecraft is in the range and range rate mode. No 
scientific data is transmitted during these intervals, of approximately 
10 minutes (8 sequences). The ranging mode was scheduled frequently when 
the spacecraft was approximately 90° 


in azimuth from the moon-sun line, 
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Experimental Observations 

In Figure 1 the results are shown for 2 orbital passes in late 
July when the moon was in the interplanetary medium. In the upper 
half of the figure the interplanetary magnetic field magnitude is 
observed to be perturbed in a very characteristic manner. As the 
spacecraft approaches the lunar optical shadow there is a small 
increase from 7.5y to 8.5y in the magnetic field at 1610 UT on July 29 
followed by a decrease to 6.5y and then a much broader increase to 
approximately 8.7y in the optical shadow. Several data points are 
missing thereafter but a decreased field is readily identified followed 
by an increased field. Then the interplanetary field returns to near 
its preshadow value in both direction and magnitude. 

It is to be noted that the perturbations of the magnetic field 

magnitude in this example are clearly identified as a characteristic 

pattern of a short dipolar 4* a longer +, and finally a short dipolar 

- + anomaly with + and - representing an increase and decrease, 

respectively. The magnetic field perturbations are about a gamma in 

an average field of 7 . 5y , a perturbation of less than 207 o . Note that 

the angular deviation of the field direction measured by 8 appears 

o 

to be no more than about 20 , while there appears to be a slight 
characteristic variation in the azimuthal direction 0 from pre to 
postshadow. 

The lower half of Figure 1 displays an additional traverse 
through the lunar wake on July 30 in which the sharply defined features 
discussed above are not identified. While a rather broad decrease is 
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observed both pre and postshadow, a shadow increase is not clearly 
detected. A unique pattern of the 9 perturbation in the shadow is 
noted. The fluctuations of the interplanetary magnetic field as 
measured by the parameters 6C and 6F are different in these two passes 
but the orientation (9,0) is approximately the same. However, there 
is no evidence for the existence of a lunar associated region in which 
the magnetic field suddenly increases at a shock surface and an onset 
of increased fluctuations is observed. For the relative position of 
the spacecraft orbit, a shock wave would be expected to be observed 
between the region defined by 0ggg = 135° to 0gg E = 45°. 

In Figure 2 are represented successive orbital passes on August 1 
in which the general pattern of fluctuations or perturbations in the 
lunar wake is somewhat different than shown in Figure 1. In Figure 2a, 
it is noted that there is a general increase of the magnetic field 
throughout the perturbed wake region although there still exists a 
suggestion of the characteristic + - + - + perturbation pattern. 
Clearly, the penumbral dipolar anomalies are not as well developed as 
in Figure la. The apparent increase in magnetic fluctuations in the 
lunar wake during this interval, as measured by 6C, is somewhat unique 
and generally not a characteristic feature of the observations . 

In Figure 2b, the well developed + - dipolar perturbation in the 
solar plasma penumbra is observed on entry but the remainder of the 
lunar wake region is almost lacking of characteristic perturbations 
detected previously. No identifiable umbral increase is noted. 

In Figure 3 are presented the 2 orbital passes on August 2 in 
which the readily identifiable dipolar anomalies in the solar plasma 
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penumbras are present although the time interval and position relative 
to the optical shadow are different. In Figure 3a, the relative 
increase in the lunar shadow of the magnetic field between 0110 and 0140 
is very small. In Figure 3b, the well-developed dipolar penumbral 
anomaly is readily identified near 1230 while the perturbation to 
the interplanetary magnetic field in the solar plasma umbra is 
missing any increase generally being less than %y. Here the 
fluctuations of the magnetic field, as measured by 6C and 6F , are seen 
to be extremely low and it is inferred that this time interval represents 
a period during which the interplanetary medium was very steady and 
uniform over a three hour interval. No evidence is found for a detached 
bow shock wave, in this, or indeed in any of these data. 

In Figure 4 are presented two successive sets of results obtained 
on August 3 and 4. The data in Figure 4a display the characteristic 
patterns of + - + - + perturbations found previously. Note here that 
the perturbations are similar to those previously obtained with readily 
identifiable increases preceding the decreases in the solar plasma 
penumbra. 

A classic example of the perturbation of the interplanetary magnetic 
field by the lunar wake is shown in Figure 4b . The characteristic 
feature of a + ■ + - + anomaly is clearly observed. In addition, there 
appear to be additional + - and - + perturbations identifiable exterior 
to those noted in the preshadow data. 

Additional observations obtained on August 4 and 5 are shown in 
Figure 5 which reveal the wide variety of perturbations of the 
interplanetary magnetic field observed in the lunar wake. In Figure 5a 
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it is seen that the leading anomalies in the solar plasma penumbra 
are much broader than previously observed but that the perturbation 
postshadow is quite different in that it is observed for a much shorter 
time interval. In Figure 5b there appears to be an absence of the 
characteristic penumbral dipolar perturbations previously observed. 

In Figure 6 are shown the last two experimental sets of observations 
to be presented, obtained on August 5 and 6. The characteristic + - + - + 
pattern of perturbations is again observed in both orbital passes. 

To summarize these experimental observations we note the following: 
generally there is an increase in magnitude of the interplanetary 
magnetic field detected in the solar plasma umbral shadow of the moon 
of up to about 3y . A relative decrease of up to 4y is 
frequently noted on either side of the umbral shadow increase. These 
decreases generally occur in or near the penumbral plasma shadow of the 
moon and are bounded on the exterior by additional increases in the 
magnitude of the field. Thus the characteristic anomalies in the 
penumbral regions consist of an exterior increase and an interior 
decrease on either side of the umbral increase. The penumbral dipolar 
anomalies are not symmetrical in their magnitude or time length as 
observed by the spacecraft. This is interpreted to represent a feature 
related to the selenocentric orbit in which the preshadow penumbral 
anomalies are observed at a radial distance of approximately 2.5 
while the postshadow dipolar anomalies are observed at only 1.6 R^. 

In addition the spacecraft is moving more slowly at the greater radial 
distance and this distorts the apparent length and extent of the anomaly. 
The perturbations observed pre-optical shadow generally extend into the 
optical shadow period. 
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A unique feature of these data is that the pre and postshadow 
perturbations of the interplanetary magnetic field consist not only 
of a decrease but also an increase of the field. Moreover, on 
several occasions there is an indication that additional increases 
and decreases exist further out from the umbra. The development of 
a pattern of alternating increases and decreases is evident in the 
data obtained in the solar plasma penumbra. Usually they appear 
in paired form, i.e., always a preshadow 4- and - while postshadow 
they appear as - + . In all cases the magnitude of the anomalies 
decreases rapidly as the order of the pattern increases. In many 
instances only the innermost dipolar anomalies are readily identified. 
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Induced Electric Currents and Associated Magnetic Fields 

In order to calculate the perturbations of the interplanetary 
magnetic field In the lunar wake, it is necessary to know the electric 
currents in the plasma. In this sector these currents and the 

perturbations of the interplanetary magnetic field are discussed 
theoretically. A first order solution of Maxwell* s equations has been 
obtained which approximates a self-consistent solution for the field 
and the currents. This solution explains both the increase of 
magnetic field in the umbra and a decrease in the penumbra. It also 
suggests the presence of periodic increases and decreases outside the 
umbra. 

The variation of the magnetic field is governed by Maxwell's 
equations 


V - B = ° (2) 

and 

V x B = J. (3) 

Introducing a vector potential A such that _B = x A, and requiring 
that V*A = 0, then we have 

V 2 A = - fUL J (4) 

~~ c ~~ 

Once the distribution of current density in the lunar wake is known, 
equation (4) can be solved for A and thus B obtained. 

In this study three sources of current are considered; the 
magnetization current, the gradient drift current and the curvature 
drift current. The magnetization current is produced by the gyration 
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of the plasma particles in the magnetic field and its density, J , 
may be expressed as 


^ - V x ( B) 


(5) 


Here < p, > denotes the average magnetic moment and N the number density 
of charged particles • Due to a gradient in the magnetic field or a 
curvature of the field lines, particles of opposite charge will drift 
in opposite directions. These drifts produce the gradient drift current, 
and the curvature drift current, J^. They may be expressed as 


and 


4 ? 



< H > B x VB 



JL x Cfe • ) 


} ( 6 ) 


where 


p = r nv 2 fd 

il -j v M 


is the pressure along the lines of force and b the unit vector in the 
direction of the magnetic field. Thus the total electric current in 
the lunar wake is taken to be 

^MGR = ' J -M + + J -R 

and equation (3) becomes 


V x B = 


4tt 

c 


4rr 

c 


MGR 


( J M + J ( 


•Jr) 


(7) 


Substituting (5) into (7) yields 
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v X C(1 + <ML>) B] = itn (T + J ) (8) 

— g c jx 

Since the magnetic moment <|jl> is an adiabatic invariant, the left- 
hand side of (8) can be expanded into three terms: 


(1 + itHM <^>) V x B - ^<(i>BxVN + ^ <ii> B x VB (9) 

B ~ ~ B B 2 

The third term in (9) is identically equal to JjQ and hence shows 
that equation (8) can be reduced to 


V x B = 


4tt 


B+4ttN< tt> 


Bx[< h > V N + — >*■ b.Vb ] 


( 10 ) 


Thus the density of induced current in the lunar wake can be expressed 
as 


~’ MGR B+4ttN< p, > 


B x[<U> V N 


+ -IL 

B 


(ID 


The average magnetic moment < 'M-> can be expressed in terms of B 0 /N q as 

0 B 

<|i> = — - — Q 

8n N o 

where 0 is the ratio of the perpendicular pressure to the magnetic 
pressure in the plasma flow upstream of the moon. For representative 
values of T x = lO^-lO 5 ^, n = 2-5 p/cirt 3 , B = 4-7 gamma in the 
interplanetary medium, 0 is of the order of unity. 
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In the lunar wake, the electron density can deviate from the ion 
density only very slightly because of electrostatic forces. The 
ion flow in the wake region has been studied by 

Whang (1967) „ In his mcd el the flow of magnetized solar wind plasma 
in the vicinity of the moon is treated as a free particle flow in the 
guiding center limit. In order to write his solution of ion density 
distribution in selenocentric solar ecliptic coordinates we consider 
that the solar wind velocity is along the heliocentric radial direction 
(-Xgs E ) anc * th e interplanetary field upstream is parallel 
to the plane of the ecliptic, the X-Y plane. The flow field is 
symmetrical about the Z = 0 plane. The motion of each guiding center 
particle is restricted to a Z = constant plane. The presence of the 
moon does not perturb the flow conditions in the region j Z |> R^. 

Dividing the perturbed flow into three regions as shown in 

Figure 7, the density distribution, N, and the parallel pressure, P , can 

be expressed as 


in region I N = (N 0 /2) erfc (y T S) 

P = P LN/N + JH YjSexpC-Yj-S 2 )] 

II II o ° 

in Region II N = (N Q /2) erfc (y 2 S) (12) 

p = P [N/N + Jn Yi S exp ( ~y| S 2 )] 

II II 

in region III N = (N 0 /2) [erfc (YiS)+ erfc(YsS)] 


p = p [N/N -h/FF Yi S exp(=Y:i . 2 S 3 )-hJ~vi y z Sexp( -y z * S° ) 1 

II 11 ° o 


where 


Yi = sin (a + X) /sin(0-cv-X) 

y z = sin (cv-X)/sin (0-ky-X) 

a = arc sin[(R M s - Z 2 )^ (X 2 + Y 2 )"^] 


and 


\ = arc tan (Y/X) 

S = U/ (2k Tu/Mj)* 
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Here T j( represents the average component of the thermal velocity of 
the ions parallel to the interplanetary magnetic field. The density 
and the parallel pressure are functions of the speed ratio, S, and the 

direction angle of the undisturbed magnetic field 0. 

Since S is the ratio of the solar wind velocity to the ion thermal 
velocity parallel to the field lines it can vary considerably 
in the interplanetary medium. For solar wind velocities of 300 to 
500 km/sec, a temperature of 2 x lO^K implies speed ratios of 16.7 
to 27.7 while for the same range of solar wind velocities and a 
temperature of 2 x 10 s °K the speed ratios range from 5.4 to 12.5 
respectively. Thus typical values of S range from about 5 to 30. 

Because higher temper atures are usually found at the times of higher 
velocities a value of S = 10 is taken to be representative of normal solar 
wind conditions, based on the speed ratio of the ions only. 

Representative lunar wake models resulting from application of 
(12) are shown in Figure 8 for speed ratios 4 and 10. The interplanetary 
magnetic field is assumed to be oriented at the theoretical spiral 
angle of 0 = 135° at 1 AU. It is seen that a large speed ratio leads 
to a longer and more narrow wake region. The projection of the sample 
Explorer 35 trajectory in Figure 8 has been chosen in the middle of 
the time interval covered by the experimental data discussed in the 
previous section. Whang (1967) has shown that the surfaces of constant 
ion density and flux are very close to being spatially coincident. Thus 
for the remainder of this paper, we shall discuss the lunar wake in 
terms of the density contours since they are explicitly related to 


the induced electric currents. 
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In this figure and in the subsequent discussion the umbra of 
the solar wind plasma shadow is defined as that region behind the 
moon for which the density value is below a particular threshold 
selected. In the case of the data shown here the value N = .01 N Q 
or 1 percent of the unperturbed density has been chosen. 

With this definition the length of the umbra in the case S = 4, 

0 = 135° is seen to be 3 while in the case S = 10 , 0 = 135° the 
length is about 8 R^. 

The length of the lunar umbra depends critically not only upon 

the speed ratio S but also field orientation 0. In Figure 9 is shown 

a parametric set of curves yielding the length of the lunar plasma 

umbra* depending upon the threshold flux level chosen, as a function 

of speed ratio S and field orientation 0. This figure illustrates 

the essential features of the lunar wake dependency upon field 

orientation. The shortest umbra is obtained when the field is oriented 

o 

perpendicular to the solar wind flow (0 = 90 ) ; an infinitely long 
umbra would be obtained for field lines parallel to the flow (0 ~ 180 ) 
Intermediate lengths are obtained for 90° < 0< 180° with the length 
of the umbra being a sensitive function of both field orientation as 
well as speed ratio S. 

Considerable symmetry exists in the theoretical solution of solar 
wind flow past the moon. This can be noted by observing that the 
sense of the field line at a given direction is unimportant. Thus the 
ion density contours shown in Figure 8 are equally valid for a field 
orientation of 0 = 315° ( = 180° + 135°). In addition, the solution 

can also be considered valid for a field orientation of 0-45 (360 
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315°) or 225° (360° " 135°) if the selenocentric solar ecliptic 
coordinate system is rotated about the X axis so that Z is now 
directed southward to the ecliptic plane. Relative positions of a 
satellite traversal of the lunar wake would be correspondingly 
rotated. 

Using the induced current density, ion density and parallel 
pressure given above (Equations (11) and (12)), numerical solutions 
of equation (4) for the perturbation of magnetic field on the ecliptic 
are obtained by an iteration process for various values of p^ Q /Px 0 » S, 

0, and 3 . The numerical method employed is discussed in detail by 
Whang (1968). 

Numerical solutions for various values of p /p (the ratio of 

II o r ± o 

parallel pressure to perpendicular pressure in the unperturbed solar 

wind) will be discussed first. The measured values of p /p vary 

II o r j_o J 

from lo5 to 4 under relatively quiet conditions (Hundhausen et al. , 1967). 

Specific numerical solutions for S * 10, 0 « 135° and 3 = 1 are presented 

in Figure 10 at two radial distances (R = 2.5 R^ and R = 4.5 R^) . The 

parameter p /p measures the relative importance of the curvature 
U° -L° 

drift current, Jr/^hGR* From this figure it can be seen that the 

curvature drift does enhance the magnitude of the magnetic field perturbation, 
but the general features of the field perturbation remain unchanged. This 
indicates that the important features of the field perturbation are 
actually controlled by the magnetization current and the gradient 

current (probably the effect of J q is itself small compared to that 
of Jj^) o Letting denote the sum of and , equation (10) shows that 
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J]ylG P ro P ort ionai to B x vN. On the plane of symmetry is 

perpendicular to that plane but points in opposite directions on either 
side of the shadow. The magnetic field induced by causes the 

increase of magnetic field in the umbra and the decreases in the penumbra. 
Specific solutions for p ^/p q ~ 3, S*7 , 3 = 1 and R = 3.5 
are presented in Figure 11 for three values o f 0, 105°, 135° and 165°. 

The noticable effect associated with the varying field orientation may 
be understood as follows. In the regions where the magnitude of yN is 
relatively large, the vector VN is nearly perpendicular to the direction 
of the solar wind* When the acute angle between the field line and the 
solar wind direction decreases, J^q increases. Increasing causes 

both larger increases and larger decreases of the magnetic field in the 
wake region as shown in Figure 11. Also the length of the umbra changes 
with 0 so that at a fixed distance from the moon the relative position 
in the lunar wake depends on 0 . Thus at a distance of R = 3.5R^, the 
scale distances of the three solutions shown in Figure 11 are 3. 5/4*0 =.88, 
3. 5/7. 8 = .45 and 3.5/14, = .25. It is seen that the largest anomaly 
is obtained for the smallest scale distance or equivalently for the 
largest 0. Note that there are indications in the bottom curve of 
relative magnitude increases exterior to the penumbral decreases, as 
well as variations within the umbra. 

Numerical solutions for different speed ratios (8=4 and 10) at 
two radial distances (R = 2 R^ and R = 4 R^) are plotted in Figure 12. 

From this figure we can see that the variations of the magnitude 
perturbations are small for different speed ratios _Lf the. field 
orientation is held constant. In this figure it may be noted that the 
two decreased regions of the penumbra become narrower and closer to each 
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other as the speed ratio S increases* The reason for this is that as 
S increases, the wake becomes narrower and the plasma density changes 
over a narrower region* The effect becomes more apparent as the scale 
distance R increases* 

In the results presented in Figure 13 the effect of varying 0 on 
the perturbation of the magnetic field can be seen* The three curves 
for 6 = 0.5, 1, and 2 are nearly similar to one another, and the 
magnitude of perturbations is approximately linearly proportional to the 
p values. As 0— 0 , which corresponds to <\i> = 0, the trivial solution 
of equation (10) is that B = everywhere and the interplanetary field 
will not be perturbed by the lunar wake* Once a solution is obtained 
for the perturbation of magnetic field for 0=1, this similarity relation 
may be used to estimate the magnitude of perturbations for other 0 
values. That is, denoting by the total field magnitude for the case 
0 = 1 and by B 0 the unperturbed interplanetary value 

B = B + H0 (B.-B ). 

o 1 o 

Where n is found to be approximately 0.90. 

Qualitatively and quantitatively these theoretical calculations 

do produce some of the features observed in the satellite data. The 
following section will discuss the relation of the calculations to the 


measurements . 
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Discussion of Results 

It has been shown in the previous section that a first order 
solution of the perturbed interplanetary magnetic field in the lunar 
wake displays characteristic umbral and penumbral anomalies similar 
to those experimentally observed. The magnitude of the umbral increases 
are theoretically predicted to be on the order of 10 to 40 percent. 

These magnitude anomalies are in good agreement with experimental 
observat ions . The innermost penumbral decreases theoretically 
predicted are less by a factor of 2 to 4 but experimental observations 
indicate they are of approximately the same magnitude as the umbral 
increases. Finally the experimental evidence reveals the existence 
of exterior penumbral increases whose effects are about the same 
as the umbral increase and the penumbral decreases. Only in the 
extreme case of the field orientation approximately parallel to the 
solar wind velocity, as shown in the bottom of Figure 11, has a small 
penumbral increase been predicted. It is believed that the failure 
of the present model to adequately explain the penumbral dipolar anomalies 
is related to its approximate nature. 

The solutions which are obtained by Whang (1968) are approximate 
and valid only in the symmetry plane of the magnetized plasma flow 
around the moon. As previously noted the induced electric currents 
are not axially symmetric . Thus if the satellite passes 
through the lunar wake at various angles with respect to the plane 
containing the magnetic field and the solar wind velocity vector, the 
magnetic perturbations observed would be expected to be different. 

Indeed the first order perturbation theory based upon the zeroth 
order model of the lunar wake suggests that the narrowest penumbral 
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anomalies are expected to be observed as the satellite passes through 
a plane perpendicular to this symmetry plane. The broadest anomalies 
would be observed in the symmetry plane itself with intermediate 
orientations and positions yielding intermediate cases between these 
two extremes • 

An attempt has been made in the presentation of the experimental 
data to consider this fact. In the upper right hand corner of each 
of the satellite traverses shown in Figures 1-6 is the satellite 
orbital trace projected into the Y-Z plane. The Z axis is selected 
perpendicular to the symmetry plane which is defined 
in terms of the average preshadow interplanetary magnetic field 
orientation and the moon-sun line. Without direct measurements of 
the interplanetary plasma flow velocity this symmetry plane is only 
approximate but should be valid within approximately 10°. In addition 
the interplanetary magnetic field is observed to fluctuate over a 
similar range and hence the diagrams only indicate approximate relative 
positions. The arrow plotted represents the trace of the satellite 
orbit through the optical shadow with the tail corresponding to the 
point of entry into the shadow and the head corresponding to the point 
of exit from the shadow. 

As can be noted in Figures 1-6 only a very few cases satisfy the 
requirement of the satellite trajectory being approximately in the 
symmetry plane of plasma flow around the moon. However, in several 
of the cases shown the field is nearly parallel to the moon-sun line 
so that the wake is nearly axially symmetric. For those cases near the 
plane of symmetry or with axial symmetry it is possible to directly 
compare the results with theory. Certain parameters are still unknown 
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however, namely the speed ratio, S, and the value of 3 . In a 
further attempt to interpret the experimental measurements the 
appropriate value of the planetary magnetic activity index Kp is 
included in Figures 1 through 6, It is known that high values of 
Kp generally correspond to increased values of solar plasma velocity, 
implying a large speed ratio, S. Unfortunately, the relationship 
of the solar plasma temperature to the activity index Kp is not yet 
known. However, unless the temperature increases much more rapidly with 
Kp than does U Q , it is still reasonable to assume that in fact S itself 
will increase because of the square root dependency of S on • 

A review of the experimental results in the light of the 
discussion of the satellite motion relative to the symmetry plane 
and the magnetic activity index leads to the following conclusion. As 
the Kp index increases the width of the wake region through which the 
satellite passes appears to decrease if due regard is given to the 
satellite motion relative to the symmetry plane. More precise 
comparisons of the theoretical model presented herein and the observed 
interplanetary magnetic field perturbations will require the introduction 
of specific plasma parameters defining the speed ratio and the 3 value 
for the plasma. 

In the work of Colburn et al. (1967), it has been suggested 
that the increases and decreases of magnetic field in the lunar wake can 
be explained by means of the equation 



4- n k T x = 


constant 


(13) 
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This describes the balance of total plasma and field pressure in a 
direction transverse to the magnetic lines of force. In general the 
magnetic field is not parallel to the direction of the wake, the magnetic 
field stress is a tensor quantity, and its magnitude is quite different 
from B 2 /8tt in a direction other than transverse to the lines of force. 
Specifically the magnetic field exerts a tension of magnitude B 2 /4 tt 
directed along the lines of force as well as a pressure. 

The argument by Colburn et al. (1967) that the penumbral decrease 
of B is due to an increase of Tj. implies that the magnetic moment is not 
conserved in the wake region. The use of an expansion fan model of the 
penumbral decreases in the work by Johnson and Midgley (1967) and Miche-1 
(1967b) is based upon the assumption that the solar plasma behavior in 

the wake is described by a fluid flow model. 

In the work presented in this paper, Maxwell's equations are solved 
directly so that the tensile stress is not ignored. It is suggested here 
that the induced electric currents due to diamagnetism of the plasma, 
the gradient drift and the curvature drift are the primary mechanisms 
disturbing the interplanetary magnetic field in the wake region. Furthermore, 
since the model is based on the independent motion of particles in their 
self consistent fields, collective effects such as a tail shock predicted 
by Michel ( 1967a, b) are not expected. Thus far the experimental results 
give no evidence for a tail shock. 

When comparing the GSFC magnetometer data with that published 
by Colburn et al. (1967) from the ARC magnetometer experiment 
the sharp decrease to 5.3y shown in their Figure 3 near 1425 UT 
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on July 31 is not observed - (See Figure 14) . From 1413-1430, 

the NASA-GSFC experiment measures a 6.9y magnitude average, which 

is approximately 0.5y less than the ARC value while the weakest field 

measured is 6.5y. It would appear that the ARC data for this 

particular interval include some spurious data points not eliminated 

in their analysis. The associated physical interpretations are inapplicable 

One particular advantage of the physical model proposed here is 
its ability to numerically predict the penumbral increases and decreases 
experimentally observed so clearly. It is suggested that the penumbral 
variations be viewed as periodic structures which arise from the 
gradient of plasma density in the disturbed flow field around the 
moon. The first order solutions produce primarily decreases in the 
penumbra but suggest a periodic structure in some cases. Further 
iteration may be expected to produce both increases and decreases 
as the self-consistent solution is approached. Work is presently 
in progress along these lines and will be reported in the future. 

In this study, the explicit effect of the induced electric fields in 
the wake region has been assumed to be negligible. If this is not valid, 
however, there will not be major differences in the solar wind 
flow pattern or magnetic field perturbations. The principal effect 
would be for the electrons, which have a much lower speed 
ratio, S e , to attempt to fill in the lunar umbra faster than the ions. 

The resultant electric field would force the ions to fill in the umbra 
faster than predicted based upon their own speed ratio, S^. Thus, the 
net result is to require the computation of an effective speed ratio, 

S p « S* < S . 
c x 


S* 
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which is lower than S^, for determining the geometry of the lunar 
wake. In the absence of direct measurements of the solar wind 
electron temperature, T e , this would then indicate a way of estimating 
T e from the discrepancy between the observed wake geometry and the 
wake predicted using a speed ratio based only upon ion temperatures. 
However, as noted in the results presented in Figure 12, this may not 
be a sensitive method to measure T e since the perturbed interplanetary 
magnetic field differs only slightly from a value of S=10 to S=4. 

An important feature of the observed magnetic perturbations in the lunar 
wake is their variability. Sometimes they are quite prominent while at 
other times they are too small to identify. For instance in Figure la 
the increases and decreases are readily identified while in lb they are 
quite vague. Similarly the passes in Figure 2, particularly 2a, show 
very small amplitude changes with only a hint of the characteristic 
+ - + - + behavior. 

These changes in the amplitude of the magnetic perturbations are 
related to the parameter 3 • As 3^ the amplitude of the perturbations 
also vanishes. Thus it is expected that the period from July 30 to 
August 1 was one of small (3 and that both before and after this time the 
plasma energy density was greater. Future comparison with the data from 
the MIT plasma probe should verify this point. Ultimately it is planned 
that measurements of both S and 3 will be used to compute the magnetic 
perturbations for the purpose of comparing with the magnetometer measurements. 
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Summary 

Repeated measurements from selenocentr ic orbit by the Lunar Explorer 
35 NASA-GSFC magnetic field experiment give no evidence for the existence 
of a lunar bow shock wave or a standing shock behind the moon. The 
interplanetary magnetic field appears to be convectively transported through 
the moon with little distortion and hence the theoretical model proposed 
by Gold (1966) and interpreted by the Luna 10 measurements (Dolginov et al., 
1967) is not supported. Variations of the magnitude of the magnetic field 
show solar plasma umbral increases of approximately 10 to 30 percent with 
periodic variations of similar magnitude in the solar plasma penumbra. 
Variations of the direction of the magnetic field are generally less than 
2(f - 30°. A characteristic pattern of the magnitude anomalies is + " + “ + 
with the interior + representing the umbral increase and the exterior + - 
and -+ representing the periodic (here dipolar) variations in the penumbra. 

From the joint experimental and theoretical results obtained thus far, 
it appears that the magnetic field perturbations in the lunar wake are of 
complex geometry depending critically upon the solar wind velocity, density and 
temperature as well as the field orientation, 0. The speed ratio, S, mea- 
sures the relative magnitudes of the solar wind velocity and thermal ion 
velocity parallel to the field line and is a critical dimensionless para- 
meter employed in the analysis. The parameter 3 is also critical. 

It is tempting to conclude that higher order iterations of the plasma 
and field equation will lead to a more self-consistent solution and thus 
will provide better agreement with observations. At present the absence 
of definitive plasma data precludes exact comparisons of the experimental 
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results for those few existing cases in which the satellite orbit is 
approximately in the symmetry plane of plasma flow. The changing nature 
of the interplanetary medium on both short and long time scales precludes 
a statistical study. As additional experimental data are obtained and in 
particular, from simultaneous observations of the interplanetary medium, 
by separate satellites a more comprehensive comparison with theory will be 
possible . 

It is suggested that the periodic nature of the field oscillations 
in the lunar wake are characteristic of the disturbed plasma flow around 
the moon. On the basis of these results and the theoretical model pre- 
sented, it is not expected that a tail shock will be observed in the lunar 
wake, as proposed by Michel (1967a, b) 
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Figure 1 


Figure 2 


Figure 3 


Figure 4 


Figure 5 


FIGURE CAPTIONS 

Magnetic field measurements, presented in selenocentric solar 

\ 

ecliptic coordinates, obtained on July 29 and 30, 1967, by 
Explorer 35. Perturbations of the interplanetary magnetic 
field identified as relative increases (+) and relative 
decreases (-) are indicated above the magnitude results. 

The Kp index of terrestrial magnetic activity for 29 July was 
3+ 4- 3- 1+0+2- 2 4, and for 30 July was 4+ 5 5-1+1 0+ 1- 0+ 
See text. 

Magnetic field measurements, presented in selenocentric solar 
ecliptic coordinates obtained on August 1, 1967 by Explorer 
35. The Kp index of terrestrial magnetic activity for 1 August 
was 1- 1- 0+ 2- 1 2- 1+ 0+ . 

Magnetic field measurements, presented in selenocentric solar 
ecliptic coordinates obtained on August 2, 1967, by Explorer 
35. The Kp index of terrestrial magnetic activity for 2 August 
was 0+ 0+ 0 1+ 0+ 0+ 0+ 1+. 

Magnetic field measurements, presented in selenocentric solar 
ecliptic coordinates, obtained on August 3 and 4, 1967, by 
Explorer 35. The Kp index of terrestrial magnetic activity 
for 3 August was 00+0+0+1112 and on 4 August was 
1 1- 2+ 1 3 3 2+ 2-. 

Magnetic field measurements, presented in selenocentric solar 
ecliptic coordinates, obtained on August 4 and 5, 1967 by 
Explorer 35. The K p index of terrestrial magnetic activity 
for 5 August was 3+ 1+ 1+ 2 2+ 2+ 2+ 2- 1. See Figure 4 
caption for K p values on 4 August. 
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Figure 6 


Figure 7 


Figure 8 


Figure 9 


Figure 10 


Magnetic field measurements, presented in selenocentr ic solar 
ecliptic coordinates, obtained on August 5 and 6, 1967, by 
Explorer 35. The Kp index of terrestrial magnetic activity on 
August 6 was 1-21+211-2+ 2+. See Figure 5 caption for 
Kp values on August 5. 

Representation of characteristic regions of perturbed plasma 
flow surrounding moon (Whang, 1967). 

Projection of orbit of Explorer 35 on ecliptic plane in lunar 

wake. The interplanetary magnetic field is assumed to lie in 

o 

the ecliptic plane at an angle of 135 to the earth-sun-line. 

Contour lines of uniform density at intervals of N = 0.01, 0.1, 

0.5, 0.9 and 0,99 of the undisturbed flow field N 0 are presented. 

U represents the umbral and P the penumbral regions in the flow 
field assuming a threshold of N = 0.01 N q and 0.99 N c . 

The length of the lunar wake, defined by the fraction of free 
stream density N/N Q , as a function of the magnetic field orientation, 
0, Two families of curves are shown corresponding to the two 
speed ratios S = 4, 10. 

Theoretical perturbed interplanetary magnetic field, F, relative 
to free stream value F q in the symmetry plane at R/Rm * 2.5 and 
4.5 for S = 10, 0 * 135° and 0=1. Different solutions for 
p /p lo = 2 and 5 are compared with the solution for “ °* 
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Figure 11 


Figure 12 


Figure 13 


Figure 14 


Theoretical perturbed interplanetary magnetic field, F, relative 
to free stream value F q in the symmetry plane at R/R^ = 3.5 for 
S = 7, 0=1, and three values of 0. 

Theoretical perturbed interplanetary magnetic field F, relative 
to free stream value F 0 and direction 0 in symmetry plane for 
0 = 135°, 0=1 and R/R^ = 2 and 4, Different solutions for 
speed ratios of 4 and 10 are indicated. The symmetry of these 
solutions is indicated. 

Theoretical perturbed interplanetary magnetic field F, relative 

to free stream value F q in the symmetry plane at R/Rm = 1.7 and 

4.0 for S = 7, 0 = 135° and p /p =2. Different solutions for 

11 ° 1 ° 

0 =0.5, 1 and 2 are indicated. 

Magnetic field measurements, presented in selenocentr ic solar 
ecliptic coordinates, obtained on July 31, 1967 by the NASA-GSFC 
experiment on Explorer 35. The Kp index of terrestrial magnetic 
activity for July 31 was 0+ 0+ 0+ 0 0+ 0+ 0+ 0+. 
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